Inelastic collisions of ultra-cold heteronuclear molecules in an optical trap 
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Ultra-cold RbCs molecules in high- lying vibrational levels of the a 3 E + ground electronic state are 
confined in an optical trap. Inelastic collision rates of these molecules with both Rb and Cs atoms 
are determined for individual vibrational levels, across an order of magnitude of binding energies. 
A simple model for the collision process is shown to accurately reproduce the observed scattering 
rates. 



The electric dipole-dipole interaction provides a long- 
range, tunable anisotropic interaction between polar 
molecules. This is fundamentally different from most in- 
teractions studied between ultra-cold atoms, which are 
typically isotropic and comparatively short-ranged. Fea- 
tures of the dipole-dipole interaction can lead to many 
novel and exciting phenomena, such as field-linked states 
1|, long-range topological order [2|, quantum chemistry 
3, Jj], and the possibility for quantum computation 
3, [a]. Furthermore, the presence of closely spaced in- 
ternal levels of the molecules, e.g. fi-doublet, rotational, 
and vibrational levels, presents a host of new possibil- 
ities for precision measurement of fundamental physics 
@, H, H, O, El- Producing ultra- cold samples of polar 
molecules will facilitate trapping and, thus, the required 
high densities and long observation times for observing 
these phenomena. 

Techniques such as Stark deceleration [12] and buffer 
gas cooling [l3| are capable of producing cold samples 
from a wide range of molecular species; however, the 
temperatures and densities currently attainable via these 
"direct cooling" methods are not sufficient for observ- 
ing many of the interesting phenomena mentioned above. 
Conversely, the association of ultra-cold atoms, either via 
a Feshbach [hj or optical resonance [HI , restricts experi- 
ments to a limited class of molecules - namely, those com- 
posed of laser cooled atoms. Nonetheless, these methods 
are approaching temperatures and densities appropriate 
for observing the aforementioned phenomena. 

In this Letter, we report the optical confinement 
of ultra-cold, vibrationally excited RbCs molecules in 
the a 3 E + ground electronic state, produced via photo- 
association (PA) of laser-cooled 85 Rb and 133 Cs atoms. 
We utilize the long observation times afforded by the 
optical trap to determine the inelastic scattering rate 
for specific vibrational levels of these molecules, with 
both 85 Rb and 133 Cs atoms, across an order of magni- 
tude of binding energies. We show that a simple model 
for the collision process accurately reproduces the ob- 
served scattering rates. We also extend this model to 
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FIG. 1: (Color Online) Schematic of the experiment showing 
the overlap of the DSM with the relevant beams and rela- 
tive position of the ion detector utilized in the state-selective 
REMPI detection. 



estimate molecule-molecule inelastic scattering rates and 
discuss implications for producing trapped samples of 
X 1 E(u = 0) RbCs molecules. 

The apparatus used in this work is shown in Fig. [1] 
Briefly, 85 Rb and 133 Cs atoms are cooled and collected 
in a du al-sp ecies, forced dark-spot magneto-optical trap 
(DSM) [lallfi- Using absorption imaging along two or- 
thogonal directions we co-locate the species and measure 
the atomic density, n, and atom number, N, as nju, = 
4(2) xlO 11 cm" 3 , N Rb = 9(1) xlO 7 , and n Cs = 5(1) xlO 11 
cm -3 , N Cs = 2(1) xlO 8 . The temperature, T, of each 
species in the DSM was measured by time-of-flight ex- 
pansion to be T Rb = (80 ± 25) and T Cs = (105 
± 40) /iK. Our optical trap is a quasi-electrostatic trap 
(QUEST) realized in a 1-D lattice configuration. The 
QUEST, represented by the green beam in Fig. [IJ is 
formed by focusing and retro-reflecting the beam from a 
vertically-aligned 100 W CO2 laser, operating with a 10.6 
fim wavelength. An acousto-optical modulator placed in 
the beam path allows for the rapid turn-off of the QUEST 
(r < 1 /is) and serves as an optical isolator for light re- 
flected back into the laser. The e~ 2 intensity beam waist 
at the focus is ~75 /xm, yielding a peak intensity of ^3 
MW/cm 2 and trap depths of ~4 mK, ~6 mK, and >9 
mK for Rb, Cs, and a 3 E + RbCs, respectively. Because 
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FIG. 2: (Color Online) Formation and detection processes 
for ultracold RbCs. (a) The PA process excites colliding 
atom pairs into bound RbCs* molecules, which (b) decay 
into a range of vibrational states of the a 3 £ potential, (c) 
Metastable a 3 E(«) molecules are excited to level i, then (d) 
ionized and subsequently detected via time-of-fiight mass 
spectrometry. The application of Stokes light (e) instead of 
the ionizing pulse (d) can be used to produce molecules in the 
absolute ground state X 1 E(n = 0). 

we utilize a lattice configuration for our QUEST, we are 
not restricted to trapping at the focus. We find it ad- 
vantageous to trap atoms and molecules away from the 
focus, where the trap volume is much larger. By mov- 
ing the point of overlap between the QUEST and DSM, 
we trap ^9 mm away from the focus, where the waist is 
~400 fim and the trap depths are reduced by ^30. In ad- 
dition to providing a larger trapping volume, this method 
mitigates the effects of QUEST- induced light-shifts [l7| • 
The energy level scheme relevant to the RbCs forma- 
tion is shown in Fig. The PA laser has an intensity of 
^2 kW/cm 2 , and its frequency is locked to an SI = 0~, Jp 
= 1+ level, located 38.02 cm" 1 below the Rb 5S 1/2 (F= 
2) + Cs 6Pi/ 2 (F = 3) atomic asymptote [l5j . Sponta- 
neous decay of this state primarily produces molecules in 
the a 3 S + state vibrational level with binding energy Eb 
— - 5.0 ± 0.6 cm -1 , to which we assign vibrational num- 
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FIG. 3: (Color Online) Typical molecular lifetime data. Here 
the number of molecules in the a 3 E + (u = vo) state with bind- 
ing energy Eb = - 5.0 ± 0.6 cm -1 is observed in the QUEST 
as a function of time. The presence of inelastic collisions be- 
tween the atoms and molecules is evidenced by the dramatic 
reduction of the molecular lifetime when atoms are present. 
With no atoms we observe molecule lifetimes consistent with 
the background gas limited lifetime seen for isolated atomic 
clouds in the trap. 



ber vq. In our previous work [18| we tentatively assigned 
t>o = 37, but this value has an uncertainty of several 
units since the depth of the a 3 E + state potential well is 
not accurately known. 

Experimental data is taken by loading the DSM for 5 s 
from background alkali vapor, provided by heated get- 
ters, in the presence of both the PA laser and QUEST. 
Atoms are not efficiently loaded into the lattice directly 
from the DSM due to the low trap depth (~100-200 /xK) 
away from the focus. Because the QUEST is substan- 
tially deeper for the molecules than the atoms (>300 /xK), 
the PA process efficiently loads RbCs into the lattice. 
Thus, molecules are loaded into the QUEST simply by 
applying the PA beam during the DSM loading process. 
From the measured atomic densities and known PA rates 
(TH . we estimate that we trap Nr^Cs ~ 10 4 molecules at 
a density of riRbCs ~ 10 9 cm -3 and temperature of T^bCs 
w 100 /iK with roug hly 7% in the a 3 £(u = v ) state. To 
study atom-molecule collisions we load atoms into the 
lattice more efficiently via an optical molasses cooling 
stage for the desired atomic species, after the DSM is 
loaded. The optical molasses stage is performed by shift- 
ing the detunings, A, of the DSM trap lasers to Am, = 
-6r, Acs — -16r for 10 ms, where T is the transition 
natural linewidth. While the hyperfine depumping beam 
of the DSM remains on during the entire molasses stage, 
the DSM hyperfine re-pumping beam is extinguished for 
the last 100 /xs to ensure that all trapped atoms are in the 
lowest (dark) hyperfine state. Loading the lattice in this 
way leads to typical densities of n^b = 2(1) xlO 11 cm -3 
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and ncs = 6(1) xlO 11 cm 3 , and temperatures of Tfy, = 
(20 ± 11) fiK and T Cs = (20 ± 15) in a volume of 
« 7r(89 /im) 2 x 960 fim, which is roughly a factor of 40 
increase in density and factor of 5 reduction in temper- 
ature compared to loading directly from the DSM. Af- 
ter the molasses stage, we apply resonant 'push'-beams 
for 10 ms to remove any undesired atoms from the lat- 
tice. After the 'push'-beam sequence, all beams except 
the QUEST are shuttered, and the molecules and any 
deliberately trapped atoms arc held in the lattice. 

After a variable delay time, the QUEST is switched 
off and the trapped molecules are state-selectively ion- 
ized using Resonance-Enhanced Multi-Photon Ionization 
(REMPI), as shown in Fig. [TJ The resulting ions are de- 
tected using time-of-flight mass spectrometry In 
this manner, we use the observed trap-lifetime of the 
molecules in the QUEST as a direct measurement of the 
molecular collision rates. 

Typical lifetime data is shown in Fig. [3]for molecules in 
the a. 3 T, + (v = Vq) state. As can be seen, the presence of 
atoms in the lattice significantly shortens the lifetime of 
the trapped molecules, which is otherwise limited by col- 
lisions with background gas. We attribute this behavior 
to inelastic collisions between the atoms and molecules. 
These losses are likely due to ro- vibrational quenching or 
hyperfine changing collisions. Each of these degrees of 
freedom carries sufficient energy that its relaxation cre- 
ates enough kinetic energy to remove both the molecule 
and the atom from the trap, e.g. one vibrational quan- 
tum is ~2 K. The number of trapped molecules, NubCsi 
evolves in time according to: 



dt 



-^atomNRbCs - ^ N RbCs- 



(1) 



Here Tbg is the loss rate due to collisions with back- 
ground gas, T a tom is the loss rate due to inelastic col- 
lisions with atoms, f3 is the molecular 2-body loss rate, 
and V is the trap volume occupied by the molecules. 
Since two-body processes are negligible compared with 
background gas collisions {PnRbCs/^BG *C 1), we use a 
fit of the data to the form of NnbCs(t) — N e~ t / T with 
' r atom + r_BGj to extract the value of r atom- r atom 

is related to the energy-dependent cross-section, a(E), 
and the relative velocity, v, as T a tom — riatom (c(E)v) , 
where ( ) denotes thermal averaging. Hence, knowledge 
of the densities and temperatures allows the determina- 
tion of the scattering rate constant, K(T) — (a(E)v), for 
these collisions. 

Because the initial PA process populates several vi- 
brational states in the a 3 E + state, we utilize the state- 
selectivity provided by the REMPI detection to measure 
data similar to that in Fig. [3] for a range of vibrational 
states and therefore, a range of binding energies, Eb ~ 



-0.5 cm to -7 cm 



The results of the collision mea- 
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FIG. 4: (Color Online) . Molecular trap-loss scattering rate 
K constant vs. binding energy Eb, for molecules in specific 
vibrational levels of the a 3 E + state. Vibrational state label 
appears below each data point. Despite more than an order of 
magnitude difference in binding energy, the scattering rates 
appear to be identical within experimental precision. The 
error bars on each point are the results of the uncertainties 
in the density and lifetime measurements. Note that the Rb 
data's error bars have been angled for clarity. The black cross- 
hatched box and red hatched box are the prediction of the 
simple model described in the text for collisions with Cs and 
Rb, respectively. The width of the boxes is due to uncertainty 
in the collision temperature. 



sured rate constants are identical within experimental 
precision, despite over an order of magnitude of varia- 
tion in Eb- Since the molecule size and vibrational en- 
ergy spacing both change substantially over this range 
of Eb, the lack of dependence of the scattering rate on 
molecular vibration hints at a unitarity limited process, 
where the details of the short-range interaction potential 
are unimportant. This view is supported by the agree- 
ment of the data with the results of a simple model of 
the collision process, shown in Fig. 2] as hatched boxes. 
This model [2Cj], which is detailed below, simply assumes 
that any collision which penetrates to short-range results 
in an inelastic trap loss event. Thus, it represents a more 
accurate estimate of the upper-bound of an inelastic pro- 
cess than the usual Langevin (unitarity-limited) cross- 
section 21], which yields a scattering rate constant that 
is roughly two to three times larger [2jJ. 

For two colliding particles, the energy-dependent scat- 
tering cross-section for the £th partial wave, with pro- 
jection m, from state i to state / in all outgoing partial 
waves £',m is generally written as 



&£,m(E, i 



t',m' 



Tl y m,e',m'{E,i — > /)| 



(2) 



surements are summarized in Fig. U as a function of Eb 
(measured relative to the a 3 E + asymptote). The mea- 



Here Tf m # m >(E,i —> f) is the 'T-matrix', whose el- 
ements represent the probability amplitude for a tran- 
sition from the incoming spherical wave, f^ im , to the 



4 



1 .50 1:40 



TABLE I: Calculated C 6 coefficients for a 3 E+ RbCs(n) 
liding with various partners given in atomic units. 



col- 



Collision Type 


(vo - 1) 


M 


(vo + 2) 


(vo + 6) 


RbCs(v) + RbCs(i)) 


65745 


65086 


64310 


61291 


Rb + RbCs 


16991 


16920 


16869 


15960 


Cs + RbCs 


19688 


19604 


19541 


18482 



outgoing wave, ^/,£' im ' and k = y // 2/iE/h 2 is the magni- 
tude of the wave-vector at collision energy E. Since our 
experiments are sensitive to the total cross-section for 
collisions that remove molecules from the trap, we must 
sum over all final states / that lead to trap loss, in addi- 
tion to the normal sum over £, m. If we assume that any 
collision that penetrates to short-range is inelastic with 
unit probability, we can re-write Eq. @ as 



a(E,i) 



ai >- 

fjjm 



,(£,*-> /) =J2^(2£+ 1)Pt(E, £), 



(3) 

where Pt(E,£) is simply the probability of transmis- 
sion to short range. We calculate Pt(E,£) by numeri- 
cally solving the Schrodinger equation for the potential 



V(r, £) 



2fir 2 



and assuming that any flux that 



does not reflect off the potential is transmitted to short- 
range and then completely lost to inelastic processes. 
This simple technique is applicable to any highly-inelastic 
process; it requires only the knowledge of the long-range 
part of the scattering potential, which in our case is given 
entirely by the value of the van der Waals coefficient, Cq, 
and reduced collision mass, /a. 

In general, the C§ constant for two colliding particles 
is given by the integral over imaginary frequency of the 
product of the particles' dynamic polarizabilities [22| . 
Since the dynamic polarizability of RbCs was calculated 
in Ref. [23| and the atomic values are well-known, the Cq 
constant as a function of vibrational level is straightfor- 
ward to calculate. These are shown for reference in Tab. 
Ill 

Using these values, Pt, <t, and K are calculated. 
The results are shown in Fig. [51 where the scatter- 
ing rate, K = (a(E,i)v), is plotted versus collision en- 
ergy for the three classes of collisions. Note that the 
average center-of-mass frame collision velocity is given 

where Ti (T M ) is 



a.s 



Sk b T ± 
7rmi 



7rm2 



8k b T„ 

TTfl 



the laboratory frame (center-of-mass frame) temperature 
and rrii is the mass of the colliding particles. Since the 
1 = 1 (p-wave) barrier heights lie at E/ks ~ 25 /jK, 15 
/iK, and 5 fiK for RbCs colliding with Rb, Cs, and RbCs, 
respectively, the collision rates fall to the values given 
solely by their s-wave contribution at the lowest energies 
on the graph. Interestingly, quantum reflection from the 
£ = potential is found to scale linearly with k as E — > 0, 
reproducing the Wigner threshold law [24| for low tem- 
perature inelastic scattering of a cx 1/fc 25]. Thus the 
calculated scattering rate remains finite at zero tempera- 
ture, despite the fact that the unitarity limited scattering 
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FIG. 5: (Color Online) Numerically calculated scattering rate 
constant vs. energy for the three types of collisions. The p- 
wave barrier height for each collision is represented by the 
vertical lines with the same plot-style (color). 



rate scales as v 1 . From the predicted molecule- molecule 
scattering rate, we expect an initial two-body loss rate of 
yN RbCs = Kn RbCs w 0.1 Hz. Given that T BG w 2 Hz, 
our measurement of background gas limited decay for the 
pure RbCs sample is consistent with the predicted scat- 
tering rate. 

In conclusion, we have demonstrated the trapping of 
heteronuclear RbCs molecules in vibrationally excited 
levels of the a 3 E + electronic ground state. Observations 
of the molecules in the trap have revealed strong inelas- 
tic collisions, presumably due to a combination of vibra- 
tional, rotational and hyperfine quenching. A simple, 
extendable model, which relies only on the knowledge of 
the van der Waals coefficient, accurately reproduces the 
observed rates by assuming that all short-range collisions 
are inelastic for all observed vibrational levels and both 
colliding atomic species. 

We are currently working towards transferring these 
trapped molecules into their absolute ground state, via 
the transfer scheme previously demonstrated in our lab 
[T^ |. It appears that with minimal improvements, e.g. 
implementing an adiabatic transfer [261 instead of the 
stimulated emission pumping used in [191 ], a sample of 
> 10 4 absolute ground state molecules X X E(?; = 0) at a 
temperature of 20 /iK and density of > 10 9 cm~ 3 can be 
created. From calculations of the adiabatic transfer with 
available laser powers, we estimate the transfer process 
can take place in ^100 /is; thus, we anticipate negligible 
loss of population due to the inelastic collisions studied 
in the present work. In fact, it appears these inelastic 
collisions could serve a useful purpose. If Cs atoms are 
deliberately loaded into the lattice, any molecule not in 
the absolute ground state will quickly (^100 ms) be re- 
moved from the trap by these inelastic collisions. By 
contrast, X 1 E(u = 0) molecules, which cannot undergo 
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inelastic collisions with Cs atoms, are unaffected by the 
presence of the atoms - even ground state molecules can 
have inelastic collisions with Rb atoms via the energet- 
ically allowed substitution reaction RbCs + Rb — ► Rb2 



+ Cs. After the excited state molecules have been re- 
moved, resonant 'push'-beams can eject the remaining 
atoms from the trap, leaving behind a pure sample of 
ground state molecules. 
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